CONTINUOUS ITERATION OF DYNAMICAL MAPS* 



R. Aldrovandi and L.P. Freitas 



Q\ INSTITUTO DE FISICA TEORICA 



O 

<U Rua Pamplona, 145 






(N 

O 
(N 

o 

Oh 



UNIVERSIDADE ESTADUAL PAULISTA 

Rua Pamplona, 145 

01405-900 Sao Paulo SP Brazil 



E-Mails: RA@AXP.IFT.UNESP.BR and LFREITAS@IFT.UNESP.BR. 

PACS numbers: 05.45.+b; 02.10.Eb 

Key words: dynamical maps, dynamical flow, iterate maps, chaotic flows 



* With support of CNPq, Brasilia. 



ABSTRACT 



A precise meaning is given to the notion of continuous iteration of 
a mapping. Usual discrete iterations are extended into a dynamical 

flow which is a homotopy of them all. The continuous iterate 

reveals that a dynamic map is formed by independent component 

modes evolving without interference with each other. 



1. Introduction 

There are two main approaches to describe the evolution of a 
dynamical system.' The first has its roots in classical mechanics - the 
solutions of the dynamical differential equations provide the continuous 
motion of the representative point in phase space.^ The second takes a quite 
different point of view: it models evolution by the successive iterations of a 
well-chosen map, so that the state is known after each step, as if the "time" 
parameter of the system were only defined at discrete values.^ It is possible 
to go from the first kind of description to the second through the snapshots 
leading to a Poincare map. Our aim here is to present the first steps into a 
converse procedure, going from a discrete to a continuous description 
while preserving the idea of iteration. This is possible if we are able to 
"interpolate" between the discrete values in such a way that the notion of 
iteration keeps its meaning in the intervals. 

Iteration is a particular case of function composition: given a basic 
map f(x) = f<i>(x) , its first iterate is f<2>(x) = [f « f](x) = f [f(x)], its n-th 
iterate is f<">(x) = f [f<""l>(x)] = f<""l>[f(x)], etc. The question is 
whether or not, given the set of functions f^'^^(x), an interpolation f^'^^(x) 
with real values of t can be found which represents the one-parameter 
continuous group (or semigroup) describing the dynamical flow of the 
system. In order to do it, f^*^-^ should satisfy the conditions 

f<t>[f<t'>(x)] = f<t'>[f<t>(x)] = f<t+t'>(x); (1.1) 

f<0>(x) = Id(x) = X. (1.2) 

We shall find in what follows a map interpolation f^'^^(x) with these 
properties. It is a well-known fact that Taylor series f(x), g(x) satisfying 



the conditions f(0) = 0, g(0) = and with nonvanishing coefficients of x 
are invertible and constitute a group by composition."^ The neutral element 
is the identity function e(x) = Id(x) = x and two functions f, g are inverse 
to each other if it is true that f [g(x)] = g [f(x)] = Id(x) = x. In this case, g = 
f^ and f = g^ (we are taking the liberty of using the word "function" 
even for purely formal series and multi- valued maps). The clue to the 
question lies in the formalism of Bell polynomials, which attributes to 
every such a function f a matrix B[f], whose inverse represents the inverse 
function and such that the matrix product represents the composition 
operation. In other words, these matrices provide a linear representation of 
the group formed by the functions with the operation of composition. 
Composition is thus represented by matrix product and, consequently, 
iterations are represented by matrix powers. Furthermore, the 
representation is faithful, and the function f is completely determined by 
B[f]. Now, in the matrix group there does exist a clear interpolation of 
discrete powers by real powers and the inverse way, going from matrices 
to functions, yields a map interpolation with the desired properties. 

Section 2 is a short presentation of Bell polynomials, with only the 
minimum information necessary to our present objective. It is shown how a 
matrix B[f] can be found which represents each formal series f, and that the 
composition f o g of two functions is represented by the (right-)product of 
the respective matrices: B[fog] = B[g]B[f]. The identity matrix corresponds 
to the identity function, B[Id] = I, the matrix B[f^'^-^] corresponding to the 
n-th iterate f^'^'^ is the n-th power B'^[f] and the Lagrange inverse f^"^^ to a 
series f is represented by the respective inverse matrix, B[f^"^^] = B"^[f]. 

The necessity of finding Br[f] for non-integer "t" leads to the problem of 



defining functions of matrices, succinctly discussed in section 3. Given a 
matrix B, there exists a very convenient basis of projectors in terms of 
which any function of B is defined in a simple way. A method is given to 
obtain the members of this basis in a closed form, in terms of powers of B. 

The procedure is applied to Bell matrices in section 4 to obtain B^ [f] = 
g[^^t>j ^Qj. ^^y y^i^g of t, from which the function f^'^^(x) can be extracted 
and shown to satisfy conditions (1.1-2). It turns out that, though it is quite 
natural to call "time" the continuous label t, this "time" is related to a 
certain class of flows, amongst all those leading to a specific Poincare map. 
There is an extra bonus: the matrix decomposition in terms of projectors is 
reflected in a decomposition of the original map, and of its iterate, into a 
sum in terms of certain "elementary functions", each one with an 
independent and well-defined time evolution. 

2. Bell matrices 

Given a formal series with vanishing constant term, 

oo 

g(x) = I ?i xJ, (2.1) 

j=l 

its Bell polynomials Bjj^[g] are certain polynomials^ in the Taylor 
coefficients g^, defined by 



Bnk(gi. g2. • • • . gn-k.i) -wi^ [g(t)]k)}^_^ . (2.2) 

dt 



Their properties are in general obtained from their appearance in the 
multinomial theorem, which reads 



OO OO ]-j 

id 'jitO ^ = I ^ B„k (gp g2' • • • ' gn-k.i)- (2.3) 

j=l n=k 

Depending on the situation, one or another of the notations 

Bnktg] = Bnk(gl' g2' • • • ' gn-k+l) = B„]^{gj} , (2.4) 

is more convenient. The symbol { g: } represents the Taylor coefficient list 
of g, with gj a typical member. Some properties coming immediately from 
the multinomial theorem are the following: 

B„i [g] = gn ; (2.5) 

B„„ [g] = (gir ; (2.6) 

B„Jcg(t)]=ckB„Jg(t)]; (2.7) 

B„Jg(ct)]-cnB„Jg(t)] , (2.8) 

where c is a constant. Given two formal Taylor series 

OO J?. OO 

f(u) = I -T uJ , g(t) = I ?i tJ , (2.9) 

their composition 

F(t) = [f o g](t) = f[g(t)] = £ ^ Fn[f; g] (2.10) 

will have the Taylor coefficients Fn given by the Faa di Bruno formula, 

Fn[f;g] = ifkBnk{gj}- (2.11) 

k=l 

Other properties can be obtained from the double generating function 



°° *n n 



e^'^^^^-l^I^IuJBnjLg]. (2.12) 

n=l ■ j=l 



Series like (2.1) constitute a group under the composition operation (f o 
g)(x) = f [g(x)]. The identity series "e" such that e(x) = Id(x) = x plays the 
role of the neutral element and each series g possesses an inverse g<-i>, 
satisfying g<-i> o g = g o g<-i> = e and given by the Lagrange inversion 
formula. The simplest example of a Bell matrix and its inverse is given by 
the well-known case of the Stirling numbers: matrices formed by the first 
and second kind Stirling numbers are inverse to each other, because they 
correspond to functions inverse to each other. In effect, consider the series 

oo . . 

g(x) = In (1 + X) = I^:^xJ, (2.13) 

j=i 
whose inverse is 

oo 

f(u) = e"-l= L J- uJ. (2.14) 

j=i 

A generating function for the Stirling numbers of the first kind Sj^Ci) is 

1 °° n 

;ir(ln(l+x))''=£ ^^s„*). (2.15) 

n=k 

It follows from (2.3) that 

Bnk [In (1+x)] = Bnk(0!,-1!,2!, -3!, ...) = 

= Bnk{(-)Ha-l)!} = Sn(k). (2.16) 

For the Stirling numbers of the second kind S^^\ the generating function is 

°° n 

^(e''-l)''= I VSn«- (2.17) 



n=k 
from which 

Bnk [e"" - 1] = Bnk (1, 1, 1, . . . ,1) = Bnk {1} = S/) . (2.18) 



The inverse property of Stirling numbers is ^Sj^^k) Sj^Ci) = 5j, the same as 

k=j 

X;B„k[ln(l+x)] Bkj[e" - 1] = 5„j . (2.19) 

k=i 

The polynomials Bj^]^[g] are the entries of a (lower-)triangular 

matrix B[g], with n as row index and k as the column index. From (2.5), 
the function coefficients constitute the first column, so that actually B[g] is 
an overcomplete representative of g. From (2.6), the eigenvalues of B[g] 
are (gO^. Triangular matrices form a group, of which the set of matrices 
(B^y) constitutes a subgroup. Hereby comes the most fascinating property 
of Bell polynomials: the matrices B[g] = (Bjj]^[g]), with the operation of 
matrix product, provide a representation of the series composition group: 

B[g]B[f]=B[fog]. (2.20) 

It is in reality an anti-representation because of the inverse order, but this 
does not represent any problem. This property comes easily by using twice 

(2.3), as 



oo oo /i.\n 

^(f[g(t)]) " = ^(I jigj) " = I ^ B„, [fl = 



j=l n=k 

CXD CXD i CXD i i 

= I I J- Bjn[g] B„k ffl = I J- I Bj„[g] B„k [f] , 

n=k j=n j=k n=k 

from which 

Bjk[f(g(t))] = t Bjn[g]Bnk[f]' (2-21) 

n=k 



which is just (2.20). Associativity can be easily checked, and it is trivial to 
see that the "identity" series e(x) = x has the representative Bjj]^[e] = d^y-, so 

that B[e] = I. Series g(t) with g^ = can be attributed a matrix, but a 

singular one and, consequently, outside the group. Summing up, infinite 
Bell matrices constitute a linear representation of the group of invertible 
formal series. If we consider only the first N rows and columns, what we 
have is an approximation, but it is important to notice that the group 
properties hold at each order N. The general aspect of a Bell matrix can be 
illustrated by the case N = 5 : 



B[g] = 



A, 











\ 


§2 


gl^ 











§3 


3glg2 


gl^ 








§4 


4gig3+3g/ 


6gl^g2 


gl' 





rr 


in^r rr 1 <r> rr 1 


^rr rr 2 , 1 fVrr 2„ 


IHrr 3„ 


„ 5 



(2.22) 



\^g5 10g2g3+5glg4 15gig2 +10gl g3 lOgl g2 gl J 



The result (2.19) is the best example of the general property 

B[f] B[f<-i>] = I . (2.23) 

It is evident that, given the series f, its inverse series can be obtained from 

B[f<-i>]=B-i[f] (2.24) 

by simple matrix inversion. The inversion properties of Bell matrices have 
been used in the study of cluster expansions for real gases.'' Because 
g[^g<n>j _ B'^[g], Bell matrices convert function iteration into matrix power 
and provide a linearization of the process of iteration. 

Suppose now that we are able to obtain the matrix B\ with t an 
arbitrary real number. The continuous iteration of g(x) will then be that 
function g^^^ such that B[g^'^^] = B*^ [g]. By (2.5), its Taylor coefficients are 



10 



fixed by g^^^^n = B[g^'^-^]ni = B^[g]ni. To arrive at B^ let us make a short 
preliminary incursion into the subject of matrix functions. 

3. Matrix functions 

Suppose a function F(^) is given which can be expanded as a power 
series F(k) = Z°°_^ Ck(^ - 'ko)^ inside the convergence circle \X - XqI < r. 
Then the function F(B), whose argument is now a given NxN matrix B, is 
defined by F(B) = Z^. Ck(B - Xo)^ and has a sense whenever the eigenvalues 
of B lie within the convergence circle. Given the eigenvalues xi, X2, . . . , 
xn, the set of eigenprojectors {Zj[B] = IxjXxjl} constitutes a basis in which 
B is written 

B = f xjZj[B], (3.1) 

and the function F(B), defined as above, can also be written' as the matrix 

F(B) = £F(xj)Zj[B]. (3.2) 

Thus, for example, eB = Z^^e^J Zj and B" = ^^i^j" Zj[B]. The basis 
{Zj[B]} depends on B, but is the same for every function F. The Zj's, 
besides being projectors (that is, idempotents, Z^j = Zj), can be normalized 
so that tr (Zj) = 1 for each j and are then orthonormal by the trace, tr(ZiZj) 

= 5ij. Other properties follow easily, for example tr[F(B)] = Z^jF(xj)Zj 
and tr[B Zj] = (xj) . If B is a normal matrix diagonalized by a matrix U, 
UBU"^ = Bji^gQ^^j, then the entries of Zk are given by (Zk)^^ = U-^^j^Uj^^ (no 
summation, of course). 
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A set of N powers of B is enough to fix the projector basis. Using 
for F(B) in (3.2) the power functions B° = I, B\ B^, . . . , B^"\ we have I 
= El,Zj; B = r1,xj Zj; B^ = £« ^xj^Zj; . . . ; B^ = £« ^xj^Zj; . . . ; bN-1 = 

^^jXj^'^Zj. For k > N, the B s are no more independent. This comes from 

the Cayley-Hamilton theorem,*^ by which B satisfies its own secular 
equation 

A(x) = det [x I - B] = (x-Xi)(x-X2)(x-X3) . . . (x-Xn) = 0. 

A(B) = will give B^ in terms of lower powers of B, so that the higher 
powers of B can be computed from the lower powers. 

Inversion of the above expressions for the powers of B in terms of 
the Zj's leads to a closed form for each Zj, 

^ (B - xi)(B - X2). ■ .(B - Xj.i)(B - Xj + i)...(B - xn-i)(B - xn) .3 3. 

i^ (xj - xi)(xj - X2) ... (xj - Xj.i)(xj - Xj + i) ... (xj - XN-i)(xj - xn) ■ ^ ■ ^ 

The function F(B) is consequently given by 

F(B) - Ej { n 1]^ } F(xj) . (3.4) 

Thus, in order to obtain F(B), it is necessary to find the eigenvalues of B 
and the detailed form of its first (N-1) powers. Though for N not too large 
the Zj[B]'s can be directly computed, we shall give closed expressions for 
them. These expressions involve some symmetric functions of the 
eingenvalues. 

Let us examine the spectrum {xj^} of B in some more detail. The 

eigenvalues Xj^ will be called "letters" and indicated collectively by the 
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"alphabet" x = {x^, X2, X3, . . . , Xj^}. A monomial is a "word". It will be 

convenient to consider both the alphabet x and its "reciprocal", the alphabet 

X* = {x*^, x*2, x*3, . . . , x*j^} where each x*. = -1/Xj. Notice that taking the 

reciprocal is an involution, x** = x. A symmetric function in the variables 

Xj, X2, X3, . . . , Xj^ is any polynomial P(xp X2, X3, . . . , Xj^) which is 

invariant under all the permutations of the Xj^'s. Only one kind of them will 
be needed here, the "j-th elementary symmetric functions", o- = sum of all 
words with j distinct letters: 

Oo[x] = 1 (by convention) 
Oi[x] = Xi+ X2 + X3 + . . . + x^ ; 

Oo[^XJ — XiXa I XiXo T ... I XiX'VT I ... I X'^Xo I X'^X^ I ... I ... I X'VT iX'Nj, 

^nL'^J ~ '^l'^2'^3 • • • '^N-l'^N • 

The symmetric functions of x and x* are related by 

(-)^cJn-jM = C7n[x]Oj[x*]. (3.5) 

Their generating function is 

£aj[x]tJ =n(l+xjt)= fld-t/xV = 5^ n(x*,-t) . 

j=0 j=l j=i j=i 

N N 

SO that Y[ (x*j - 1) = cjn[x*] ^ CJj[x]tJ . We use the involution property and 

j=i j=0 

(3.5) to write the general expression 

N N N 

n(Xj-t) = On[x] £ Oj[x*]tJ =£(-)J0N.j[x]tJ .(3.6) 

j=l j=0 j=0 
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The j-th eigenvalue is absent in the numerator of expression (3.3) for Zj. 
We shall need some results involving an alphabet with one missing letter. 
Let Oji[x] be the sum of all j-products of the alphabet x, but excluding xi. 

For example, OniLx] = n^_^j x^. We put by convention Ooi = 1 and find that 

k k 

aki[x] = I (-) PXiPOk.p[x] = I (- Xi) k-Joj[x] . (3.7) 

p=0 j=0 

In the absence of the i-th letter, (3.6) becomes 

N N 

Y[ (Xj - 1) = GNiLx] £ Oji[x*]tJ . (3.8) 

j=l;j^i j=0 



The projectors (3.3) are then 

N 



-^ ^k-B k=o 



L^kj[x*]Bl 



^j[^^ = n ^ = "W ^ (3.9) 

k=0 



clearly written in the basis {I, B, B^, . . ., B^"^}. In our application to Bell 
matrices, it will be convenient to use instead the basis {B, B^, . . ., B^}. 
This is due to the fact that we shall prefer to start from the matrix B[g], 
corresponding to g, and not from the matrix I, corresponding to the 
identity function. The mappings of interest, like the logistic map for 
example, have a general "n" aspect and cannot be obtained continuously 
from the identity map: the identity map has Brower'' degree 1, while the 

dynamical maps have degree 0. The Cayley-Hamilton theorem implies Z^q 
Oj[x*] BJ = 0, from which we obtain the identity B^ as 
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N 

I = - J^Oj[x*] BJ . (3.10) 

Replacing this identity in (3.9), the projectors are recast into another form, 

N N 

Zi[B] = ^^ = ^ . (3.11) 

1 + L ^ki[x*]x> £{aki[x*]-Ok[x*]}Xik 



k=l k=l 

Using (3.7) and (3.5) we get 

i:{I(Xiy-'^aj[x*]}Bk i:{£(Xi)i-k(-yoN.j[x]}Bk 

Zi[B]= ^^y^^ = '-^ , (3.12) 

£(N-j)(Xi)'Oj[x*] I(N-j)(Xiy(-)JON.j[x] 

j=o j=o 

where we have also used 2:^=i{2:':q (Xi)JOj[x*]} = 2:^Q(N-j)(Xi)'Oj[x*]. There 

is a good immediate check: replacing B by the eigenvalues we find, as 
expected, 

Zi[xJ = 5ikI. (3.13) 

The projectors are now clearly in the basis {B, B^, . . ., B^}. Actually, 
each Zi is now just that given in the basis {I, B, B^, . . ., B^"^} multiplied 
by (B/xi): instead of (3.3), 

^ (B - xi)(B - X2). ■ .(B - Xj.i) B (B - Xj + i)...(B - xn-i)(B - xn) 
i^ (xj - xi)(xj - X2) ... (xj - Xj.i) Xj (xj - Xj + i) ... (xj - XN-i)(xj - xn) 

(3.14) 
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The Zi's and the powers B^ can be seen as components of two formal 
column "vectors". The linear conditions B" = Z^jXj"Zj are then represented 
by a matrix A = [xj"] , 



B2 



/ ^1 ^2 % 

2 2 2 
^1 ^2 ^3 



VBV V V ^2" %" 



Xn a 



X 



N 



X 



f^A 



, (3.15) 



n" y v^n J 



and what we have done has been to obtain its inverse: 



N 



Zi[B] = l^ [A-ljikB^; 



(3.16) 



k=l 



k-1 



£(XirH-)'aN.j[x] 



It _ -±1. 



[^" ]ik- N-1 



(3.17) 



I(N-j)(Xi)J(-)ioN.j[x] 
j=o 



It seems a difficult task to improve the above expressions, as it would 
mean knowing a closed analytical expression for the recurrent summation 

of the form Z'^^qU-'Oj[x*]. A closed expression for aj[x*] would be necessary 

and, even for the simple alphabet consisting of powers of a fixed letter a, 
which we shall find in the application to Bell matrices, this would be 
equivalent to solving an as yet unsolved problem in Combinatorics. In 
effect, in terms of such an alphabet {aJ} with N letters, the symmetric 



_ yN(N+l)/2 

" J=l 

partitions of j into k unequal summands, each one < N. These conditional 



function is given by Ok = Z.^^^ ' q-y^^a), where qj,k,N = number of 



16 



partition numbers have the generating function rTj^Jl+ua^ = 1 + S^j,>^ 

q: j^ js^aJu , but have no known closed expression. They are calculated, one 

by one, just in this way. '" 

Bell matrices are not normal, that is, they do not commute with their 
transposes. Normality is the condition for diagonalizability. This means 
that Bell matrices cannot be put into diagonal form by a similarity 
transformation. As it happens, this will not be a difficulty because we know 
their eigenvalues. That functions of matrices are completely determined by 
their spectra is justified on much more general grounds. Matrix algebras 
are very particular kinds of von Neumann algebras and it is a very strong 
result of the still more general theory of Banach algebras'^ that functions 
on such spaces, as long as they can be defined, are fixed by the spectra. 
Another point worth mentioning is that the infinite Bell matrices which 
constitute the true, complete representation of the group of invertible series 
will belong, as limits for N ^ oo of NxN matrices, to a hyperfinite von 
Neumann algebra. Our considerations here are purely formal from the 
mathematical point of view, as we are only discussing formal series. We 
are not concerned with the topological intricacies involved in the 
convergence problems, though they surely deserve a detailed study. By the 
way, the infinite algebra generated by Bell matrices would provide a good 
guide in the study of function algebras with the composition operation. 

4. The continuous iterate 

We are now in condition to find, given a function g, the matrix B* [g] 
and its corresponding function, the continuous iterate g^^^(x). As many 
things - the B[g] spectrum for example - will depend only on the first- 
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order Taylor coefficient, we shall put gj = a. By (2.6), the letters in the 
eigenvalue-alphabet of B[g] will be simple powers of a and the alphabet a = 
(a\ a^,...., a^) will have the reciprocal a* = (-a"\ -a"^,...., -a"^). The matrix 
A has entries (A^y) = (a^'^) and the projectors 



N 



Zi[B] = l^ A-SkB'^Lg] 

k=l 



(4.1) 



have now the coefficients 



k-l 



£ ai^-k)(-)ioN.j[a] 



L^ Jik- N r-1 

1 1 ay(-)JoN-j[a] 

r=l j=0 

We can verify easily that tr Zi[B] = 1 and BZ^ = a^Z^. A consequence of the 
latter is 

f(B)Zi = mz^ , (4.2) 

whose particular case 

B^Zi^a^t^. (4.3) 

will be helpful later on. To give an idea of their aspect, we show the 
projector matrices Zi^^^[B[g]] for N = 3: 



7 (3)_ 



^ 1 

gl(l-gl) ^ ^ 

3g/+g3(l-gi) 



A 







r 

-§2 



. 7 (3) _ 



\g^(i-g,ni+g^) J 



gi(i-gi) 



0\ 



1 



-3g2 



3g2 



\^gi (1-gi)' Sid-Si) ) 
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f 



7 (3)_ 








0^ 




3g2 -gig3(l-gi) - 3g, 



1 



^gl2(l-gj)2(l+g^) gi(l-gi) ^ j 

As they have not the form (2.22), they cannot be the Bell matrices of any 
function. They inherit, however, a good property of the Bell matrices: for 
each N, the projectors Z:*^^^ contain, in their higher rows, the projectors 

T.^^ for all k < N. The upper-left 2x2 submatrices of the Zj^^^^^'s above are 

just Z,(2). 

We can take the first-column entries of the matrix (4.1) as Taylor 
coefficients defining the functions 



R/^^(x) = I %. [Zi[B]]ri 

-i[A-\kg<k>(x)], 



(4.4) 
(4.5) 



k=l 



To each projector Z^ corresponds such an "elementary function" R/^^(x), a 
relationship reflecting in part that between the series and their Bell 
matrices. Taking the summation of all Rj'^'^^'s in (4.4) and using I = ^[liZi, 
we find 



N 



N 



I R,""(x) = I "fr E Zi[B]],i = I \ 8,1 = X. (4.6) 

i=l 1^1 i=l 1^1 

Thus, just as the projectors give a decomposition of the identity matrix, the 
elementary functions provide a decomposition of the identity function. 
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N 

Id = J^ R/N) . (4.7) 

i=l 

Applying to B[g] the general form (3.2) for the function of a 
matrix, we have 

B^ = '£ a^t Zi[B] = X; C(N)^(t) Rk , (4.8) 

i=l k=l 

where 

C%(t) = I ^'A-\^. (4.9) 

i=l 

The coefficients of the continuum iterate function 

OO J- 

g<^>(x) = I gr<^> p7 (4.10) 



r=l 

will be 

g<'> = [Rt [g] ]^i = £ ait [Zi[B]]ri = £ C(N)j^(t) Bk[g] ]^^ 

i=l k=l 

= I C%(t) gr<k>. (4.11) 

k=l 

Therefore, 

g<t>(x) = t C(N)^(t) g<k>(x) = t[t ^' A-\^]g<^>(x) . (4.12) 

k=l k=l i=l 

Time dependence is factorized in the alternative form 

N 

g<t>(x) = £ a^tRj^^N^Cx). (4. 13) 

k=l 

One of the announced properties, (1.2), follows immediately: g^^^(x) = 
Id(x). For t = 1, a suggestive decomposition of the function comes up: 
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N 

g(x) = E ak R^(N)(^). (4.14) 

k=l 

It remains to show that (4.13) does satisfy property (1.1). Notice 
first that, from the definition (4.4) of R/^^ and the multinomial theorem 

(2.3), 



R/N)[g<t>(^)].£ [Zi[B]]ri[^^ =11 ^^ Bjr[g<t>] [Zi[B]]ri = 

1^1 1^1 j>r 

OO . • CXD . 

I ""j^ I BV[g][Zi[B]]ri = E ^^ [B^Zilji • (4-15) 
j=i ,^1 j=i 

Therefore, 

N °° r N 

g<t'>(g<t>(x)) = £ ait' R.(N)[g<t>(^)] = I ""rf I ^'' [BtZilri = 

i=l 1^1 i=l 



OO XT OO OO 



= I n [B' I a"'Zi],i = I 'i, [B' B'],i= I 'i, g<'«'>,i = g<'«'>(x), 

1^1 i=l 1^1 1^1 

just the result looked for. Notice also that, using (4.3) and (4.4) in (4.15), 
we obtain 

and, consequently, 

R^(N)[g<t+t'>(^)j = ak(t+t')Rj^(N)(x) = a^t' Rk(N)[g<t>(x)j _ (4_ ^7) 
The function decomposition is preserved in time, as we can write 

k=l k=l 

All the above expressions hold for each value of N and give 
approximations to that order. Of course, exact results would only really 
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come out when N ^ oo. The exact expressions would, as long as they are 
defined, take forms like 

Rj(x)=I[A-Vg^^^«]; (4-19) 

k>l 

g(x) =£akRk(x); (4.20) 

k>l 

g<t>(x) = £ akt Rj^(x). (4.21) 

k>l 

It is tempting to conjecture that the elementary functions R^ play, in 

function space and with the operation of composition, the role of 
projectors. In fact, this is not so. Their properties mirror only partially 
those of the projectors Z^. If we calculate Rj [Rj(x)] using (4.4) and the 

multinomial theorem, we find Rj [Rj(x)] = Rj(x-'/j!). The composition gives 
rise to a change in the variable. Thus, the functions R^ fail even to satisfy 
the defining idempotent property: in general, Rj o Rj ;t Rj. 

We can introduce £ = In a and rewrite (4.13) as 

x(t) = g<t>(xo) = I eket Rj^(xo) , (4.22) 

k>l 

a decomposition of g^'^^(x) into a sum of modes, each one evolving 
independently according to 

Xk(t) = ek«x^(0) = ek« R^(^^), (4.23) 

The imaginary "frequency" ke plays the role of a "modular Lyapunov 
exponent" of the k-th mode. We have thus a "multi-hamiltonian" flow x(t) 
with one hamiltonian for each projector component: equations (4.16,17) 
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show Rj^*^^^ as a representation of the one-dimensional group engendered by 
the k-th dynamical flow 

Rk[g<'^(x)] = e^^t Rk[g<^>(x)] . (4.24) 

The function g^'^^(x) is actually a Lagrange interpolation in the 
convenient variable (g^f = a^ which coincides with the discrete iterates at 

each integer value of t and keeps the meaning of a continuum iterate in 
between. It is a homotopy of all the usual discrete iterations. As announced 
in Section 3, we have used the basis {B, B^, B^, . . . , B^}, instead of {I, B, 
B^, . . . , B^"^}, because in the latter the matrix "I" would correspond to 
the identity function Id(x) = x, which does not belong to the same 
homotopy class. 

The variable t has been given the sense of a "time". If it is really 
time, g^'^^(x) will give the n-th point of a Poincare map. Now, there are in 
principle (infinitely) many dynamical flows corresponding to a given 
Poincare map. The function g^'^^(x) as given above would correspond to a 
class of them, that of the flows with equal intervals of time between 
successive points on the Poincare section. To consider other cases, we 
recall that any strictly monotonous function of a first-given "time" is 
another "time". 

We cannot resist exhibiting a last "thing of beauty". By the very 
manner it has been obtained, the expression (4.12) for the continuum 
iterate is equivalent to the highly mnemonic determinantal equation 
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g<l>(x) a a2 
g<'>(x) a^ a^ 
g<^>(x) a^ a^ 



g<N>(x) a^ a2N 
g<t>(x) at a2t 



a^ 

„3N 



jN2 
„Nt 



= 0. 



(4.25) 



Expansion along the first column or the last row and comparison with 
(4.11) and (4.13) will give determinant expressions for C^^^]^(t) and 

Rk^^^(x). 

5. Concluding remarks 

In the passage from functions to Bell matrices, composition is 
translated into matrix product and iteration into matrix powers. Continuous 
powers of matrices have a sound meaning, which is translated back into 
continuous iteration and makes of it also a sound concept, respecting 
properties (1.1-2). The usual discrete iterations are thus extended into a 
continuous flow. The term-by-term factorization of the x- and t- 
dependencies of g^^^(x) reveals itself as a decomposition into independently 
evolving modes, one for each projector of the corresponding Bell matrix. 

Once it is established that iteration can have a continuous meaning 
and is, furthermore, a homotopy, new possibilities are made open to study, 
such as the use of the conservation of degrees like Brower's and Morse's. 
It is our contention that, with the continuum version of iterations, it will be 
possible to get a better understanding of the detailed unfolding of 
bifurcations and of the general relationship between the differential and the 
mapping approaches to chaotic dynamics. 
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